The fourth virial coefficient of a fluid of hard spheres in odd dimensions 



I. Lyberg * 

Institute for Theoretical Physics, State University of New York, Stony Brook, NY 1 1794-3840 

(Dated: February 2, 2008) 

The fourth virial coefficient is calculated exactly for a fluid of hard spheres in odd dimensions up 
to ff. 

YITP-SB-04-50 

Keywords: hard spheres, virial coefficients 

I. INTRODUCTION 

The virial expansion of a fluid of hard spheres has been studied for more than a century, but has not yet been 
solved. Clisby and McCoy [lj recently made some progress by calculating the exact value of the fourth virial coefficient 
in even dimensions 4 through 12. Previously the fourth virial coefficient had been analytically calculated in three 
dimensions by Boltzmann in 1899, and in two dimensions by Rowlinson Q in 1964. The second virial coefficient in 
three dimensions had already been calculated by van der Waals Q and the third had been calculated independently 
by Boltzmann 5] and Jager |6j . Luban and Baram 7] found functions that give the second and third virial coefficients 
in any dimension. They also found general functions for the two lower order diagrams in the fourth virial coefficient, 
but they did not find such a function for the complete star. Luban and Michels Q later obtained the four point 
complete star in 4 and 5 dimensions analytically as triple infinite sums. 

In this paper we compute exactly the general expression for the fourth virial coefficient of a fluid of hard spheres 
in odd dimensions up to 11. 

We consider the low density expansion of a fluid of D-dimensional hard spheres. The position of the ith particle 
is written as r^, and the distance between two points is written as r.y = |r^ — r ^ | . <fi{ r ij) is the pair potential of two 
points. For a fluid of hard spheres this is 



r ij) = 



if TiA < a 



if m > 



'■' (1) 



The Hamiltonian is 



N N 2 



i<j=2 i=l 

For low densities 



ff=£<^)+£f^ (2) 



p/kT = p+J2 B ^p n ( 3 ) 



n=2 



where p = N/V. 

There are several systematic ways to formalise the calculation of the virial coefficients B n . One of these methods 
is the Mayer expansion ||, [ljj. I n this expansion a function called the Mayer / function is defined: 

/(r tf ) = e -*(«*>/w - 1 (4) 

The virial coefficients are given by 

B n+l = / - / V n+1 (r 1 ,...,r n )d D r 1 ...d D r n (5) 

n+ln\VJ v Jy 
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where V n -\-i is the collection of labelled biconnected Mayer diagrams with n points. Each bond of these diagrams 
represents a function /(r^) in the integrand of JSJ. Explicitly 

B2 = -\ J /(r 12 )d D r 2 = -\ (6) 



B * = -\j v f(ru)f(ri 3 )f(r23)d D r 2 d D T 3 = ^A (?) 



B4=-^-|0-|n (8) 



For hard spheres the Mayer / function is 

From now on we let a = 1. If a hard sphere potential is considered, B 2 and B 3 can be explicitly evaluated using the 
formulae Q 

B2 = 2TW2TI) (10) 

and 

B 3 _ 4T(l + D/2) r /3 . D 



B 2 2 7rV2r((l + J D)/2)y 



sin ^ (11) 



The integral in (|ll|l is easily evaluated. Let m be any positive integer, and let u be any positive number. According 
to reference ^1], p. 159 



(2m -1)!!, 



sin 1 di = 



2 m m! 



cosm J . 2m _i v-^ (2m — l)(2m— 3)...(2m — 2fc + 1) . 2 m-2fr-i /-.^s 

— — sin 2m 1 m+> v -T7 -T7 „ N , rr — -sin 2m 2k 1 u\ (12) 

2m I ^ 2 k {m-l){m~2)...(m-k) J v ; 



and 

2 m m! 



sin 2m+1 a: 



(2m + l)(2m- 1)!! 

cos it 



sm 2 ™ u + Y , 2*+Mm- l) (m-fc) ^2,^2,-2 „ (13) 



2m +1 I (2m-l)(2m-3)...(2m-2fe-l) 



Table [|] shows the values of the second and third virial coefficients in dimensions two to eight. The evaluation of 
-B4 using the Mayer formalism is much more difficult. This calculation was done in two dimensions by Rowlinson in 
1964 and recently in dimensions 4, 6, 8, 10 and 12 by Clisby and McCoy pj. The results are shown in table ITT1 
Explicitly, the method that was used was to calculate the volume of the intersection of three spheres 



v D (ri2,ri3,r 2 3) = ~ f(ri4)f(r 2 4,)f(r 3 i)d r 4 (14) 
Jv 

as an intermediate step. It follows from (|14|l that 

f(r 13 )f{r 23 )f{r 2i )v D (r 12 ,r 13 ,r 23 )d D Y 1 d D Y 2 d D T 3 (15) 



v Jv Jv 



Rowlinson had prevoiusly calculated v 3 (r\ 2 , r 13 , r 23 ) [12j, but no one has so far calculated the three dimensional 
complete star using <|15[) • The reason is that there are elliptic integrals in the odd dimensional case that cancel in the 
even dimensional case. 



3 



TABLE I: The second and third virial coefficients 



11 




D, / d2 
£>3/£>2 


decimal expansion 


o 
Z 


7T/Z 


4/0 — V o/TT 


n 70000 


3 


27f/3 


5/8 


0.625 


4 


7T 2 /4 


4/3 - (3/2)V3/tt 


0.50634... 


5 


4tt 2 /15 


53/2 7 


0.41406... 


6 


7T 3 /12 


4/3 - (9/5)V3/tt 


0.34094... 


7 


8tt 3 /105 


289/2 10 


0.28222... 


8 


tt 4 /48 


4/3 - (279/140)a/3/tt 


0.23461... 



The history of the computation of JB4 in 3 dimensions dates back to the end of the nineteenth century [13J . Van 
der Waals formulated a sum of integrals which he thought would give B4. However, there was one integral which he 
could not evaluate (This was the one which is today called the complete star). Van Laar managed to evaluate this 
integral and published his result in 1899 14]. Boltzmann contested van der Waals' formulation of the problem, and 
using the correct virial series expansion he published the correct result in the same year |2j. Boltzmann's result was 

B 4 _ 2707 219 \/2 4131 arccos(l/3) 

~Bj ~ 4480 + 2240 V ~ 4480 tt ^ ' 

This result was confirmed in 1952 by Nijboer and van Hove fl5| using what is called the two center formalism. 

We shall extend this two centre computation to odd dimensions through D = 11. The analytic results are given in 
table El 

The virial coefficient B4 in odd dimensions has previously been computed by Monte Carlo methods by Ree and 
Hoover 

and Clisby and McCoy [13 . These numerical results gave the first demonstration that the hard sphere 
virial coefficients can be negative. The question of negativity of hard sphere virial coefficients is of great theoretical 
importance, and in dimensions D < 4 Monte Carlo investigations have thus far seen only positive B n for n < 10. 
The reliability of the conclusions relies on the precision of the Monte Carlo evaluation and therefore it is of interest 
to have an independent confirmation of the error quoted with the Monte Carlo result. Hence it is most interesting 
to compare the exact result of B4 for D = 5, 7, 9, 11 with the previous Monte Carlo results. This comparison is 
made in table EJ Here we see that the exact results lie within 3 times the quoted error of reference ^7|- This gives 
some measure of the accuracy of the Monte Carlo results of reference The exact results for D = 4, 5 are within 
the errors of the computation by reference |l8[ . On the other hand, the estimated value of the infinite triple sum for 
D = 5 given by reference [j| is far outside the stated error. 

In section II we review the relation between the two center formalism and the Mayer formalism. In section III we 
use the two center formalism to evaluate B4 in dimensions 5, 7, 9 and 11. We conclude in section IV with a discussion 
of the prospects of obtaining B5 and other higher order coefficients exactly. 



II. THE TWO CENTER FORMALISM 

The two center formalism was invented by de Boer in 1949 This formalism is equivalent ^(| to the Mayer 

formalism, and in the case of hard spheres it especially useful since it allows the reduction of the dimension of the 
integral by D. The invention of this formalism is what inspired Nijboer and van Hove to confirm Boltzmann's result 
for Bi in 1952 

According to the Mayer formalism, B4 is given by ©■ According to the two center formalism 

S = -4s 2 M(i) 

□ = -fs a n(i) 

Thus 



b 4 = b 2 Q n(i) + \v la) + 2 n(i) + n(i)) 



(18) 
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TABLE II: Exact and numerical values of the fourth virial coefficient 



D_ 

2 

3 
4 



Bi/B 2 



decimal expansion 



2 - + loX 

Z TT TT* 



2707 i 219 V2 _ 4131 arccos (1/3) 
4480 ' 2240 tt 4480 tt 



27 VS , 832 1 

4 TT ' 45 TT 7 



25315393 i 8888425 _ 07183425 arccos (1/3) 

32800768 ' 16400384 tt 32800768 tt 



6 

7 

8 

9 

10 

11 

12 



81 s/3 , 38848 1 
10 7T ' 1575 TT 7 



299189248759 , 159966456685 y/2 _ 292926667005 arccos (1/3) 
J<!O-,96061184 + 435894091776 tt 96865353728 tt 



2511 V3 17605024 1 
280 7r 



9 _ 2511 V3 i 

z oan „ T (i()(i375 jr 7 



2886207717678787 , 2698457589952103 x/2 _ 865606(i770O.vi~i2:{ arccos (1/3) 
2281372811001856 '5703432027504640 tt 2281372M l"Oi vl0 tt 



O _ 2673 73 , 49048616 1 
280 I ^ 1528065 It 1 

17357449486516274011 , 16554115383300832799 V2 _ 52251492946866520923 arccos (1/3) 
11932824186709344256 ' 29832060466773360640 tt 11932824186709344256 tt 



9 _ 2187 s/3 , 11565604768 1 
Z 220 tt T 337702365 ^ 



0.53223180... 

0.28694950598... 

0.15184606235... 
0.151846054(20) [8] 
0.15184(7) [18] 
0.07597248028... 
0.075972512(4) [8] 
0.07592(6) [18] 
0.075978(4) [1JJ 
0.03336314... 

0.00986494662... 
0.009873(3) [17] 
-0.00255768... 

-0.00858079817... 
-0.008575(3) [17] 
-0.01096248... 

-0.01133719858... 
-0.011333(3) [iJJ 
-0.010670281... 



Here the circles indicate points that are not integrated over, and the number 1 indicates that the distance between 
these two points is 1. We shall use the same notation as Nijboer and van Hove [l^. Thus 



X(ri2) 
(Si(ri 2 )) 2 

<p(ri2) 



= N 



(19) 



The functions gi(ri 2 ), <f(ri2) and V'( r i2) are easily calculated for any D. The calculation of these diagrams in three 
dimensions is described in the paper by Nijboer and van Hove 0|. It is easy to do the same calculation in higher 
odd dimensions, but we shall omit this since the lower order Mayer diagrams in are already known in terms of 
hypergeometric functions. Luban and Baram showed that 



□ 



5-D+4 



r(j> + i)[r(£>/2 + i)] 3 

tt r(3£>/2 + l)[T((£> + 3)/2)]- 



1 



-D+l D + 3 D + 3 



and 



where 





B 2 



l D*[T{D/2)f / dy y[g D/2 (y)f 



9v{y) 



dx x u [J u {x)] 2 J v -i(xy) 



If D is odd, then according to reference [Tj, p. 1071 



3F2 



-D+l D + 3 D + 3 



n(n — l)...(n — k + 1) 



fc=o 



2 fe [{n + k + \){n + k)...{n + 2)f 



(20) 



(21) 



(22) 



(23) 
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where D = 2n + 1 and (-1)!! = (-1)° = 1. Joslin [U pointed out that 



g tyS = J r(4l/2)"r(l/2) ^arcsin (y/2) *P C0 ^ ^/ 2 ) if V < 2 (24) 

1 if y > 2 

Thus, if y < 2, n is a positive integer and D = In + 1, then 

2 -D/2 y D/2-l 

gD ' 2{y) = r(n + i)r(i/2) 2 x 

7 1 2 "" ! <») 



I 2n + 1 (2n - 1)!! 2n + lT va/ /y ^ (2n - l)(2n - 3)...(2n - 2fc - 1) 

Thus using l|23[l and (|24|l . the expressions in (|20(l and i|21|) can be explicitly computed in odd dimensions. Clearly 
both of these are rational numbers in odd dimensions. 

III. INTEGRATION OF THE COMPLETE STAR 

We aim to obtain a general expression for x(l)- The only dimensions lower than 12 for which the exact result has 
not been published before are D — 5, 7, 9, 11. We shall calculate \ in dimensions D = 2n + l. When n is an integer, 
D is an odd integer. However, n need not be an integer. If n is a half integer, then the calculation below is still valid 
and gives B4 in even dimensions. If n is some other positive real number, then the calculation below may be used to 
obtain B4 in continous dimensions. We will use the fact that 

na > 1 (26) 

According to 119(1 

X(r i2 ) =11 f(r 13 )f(r u )f(r23)f(r 2 4)f(r 34 )d D r 3 d D r 4 (27) 
Jv Jv 

We define 

F(h)= I fin^^-^dPvi (28) 
Jv 

where h = |h|. It can be shown that 7] 

F W = -^JoM^h) (29) 
where J„ is a Bessel function of order v. We define 

G(h,r 12 )=/ f(r 13 )f(r 23 )e 2 ^ h ^~^+^d D T 3 (30) 
Jv 

Clearly 

X(rm)= f F(h){G(h,r 12 )fd D h (31) 

In D dimensions, we write r = (xi,X2, ...,Xd-i, z) = (x, z) and h = (h. x ,h z ). ri and r 2 are placed on the z axis in 
such a way that I - ! + r 2 = 0. From now on, r 12 will be written as r. We first simplify G(h, r). According to l|3U|) 

G(h,r) = 2f dz cos(2n zh z ) f d 2 "x f([x 2 + (z + r/2) 2 ] 1 / 2 )e 2 " h ^ x (32) 

^0 v/V'n{r|z=constant} 

where x = |x|. The integral over the hyperplane {r\z = constant} in l|32|) has the same form as the integral in i|28|) if 
D is replaced by 2n. It therefore follows from (|29|1 that 



2 f 1 -''/ 2 

G(h,r)=- — / dzcos(2^ 2 z)[l-(r/2 + z) 2 ] n / 2 J„(2^[l-(r/2 + z) 2 ] 1 / 2 ) (33) 

11 x Jo 



where h x = |h x |. I|31[l can be rewritten as 

x (r) = dh z [ d 2n h x F(h)[G(h,r)} 2 (34) 

J-oc JR 2 » 

Since F(h) and G(h,r) are spherically symmetric in the hyperplane {h|/i z = constant}, 13411 can be simplified as 

/oo poo 
dh z dh x h^FWlGfar)] 2 (35) 
-oo JO 

where VL 2n -i = area(S 2n " 1 ) = It follows from ®, ESI and that 

o n fl-r/2 fl-r/2 

x ( r ) = _° / dz[l- (r/2 + z) 2 ]"/ 2 / <fo'[l - (r/2 + z') 2 ]" /2 



n 



r (") Jo JO 
dh*Y J n (2tt[1 - (r/2 + z) 2 ] 1/2 ^) Jr. (2tt[1 - (r/2 + z') 2 ] 1/2 ^) 

dk * (h2±h2\D/4 J D/* ( 27r ( h x + ^) 1/2 j C ° S ( 27F ^ Z ) C0S (27T^«') (36) 
oo y^x "T" ">z) 

We rewrite cos (2nh z z) cos (2nh z z r ) as 

cos (2nh z z) cos (2nh z z r ) = ^{cos (2nh z (z + z')) + cos(2irh z (z - z'))} (37) 
According to reference p. 772 

= i-[l - (z ± z') 2 ] n/2 ^ (27r^[l - (z± z') 2 ] 1/2 ) (38) 



Thus x(r) can be reduced to a three dimensional integral. So 

47T" 



4 n /-l-r/2 /-l-r/2 

x ( r ) = -__/ dz(a(r/2,z)/2Tr) n dz'(a(r/2, z')/2-k) j 

o Jo 



oo -i 

dK—^ J„(a(r/2, z)h x )J n {a(r/2, z')h x ) {(a(z, z')/2ir) n J n (a(z, z')h x ) + {a(z, -z')/2n) n J n (a{z, -z')h x )}(39) 
where 



o 



a(z,z') = 2ttV1 - (z + z') 2 (40) 
Now we have to evaluate the integral / given by 



/ / J n {a{r/2,z)x)J n {a{r/2,z')x)J n {a{z,z')x)-^ l dx (41) 



We integrate by parts and use the recursion relations for Bessel functions 



and 



Ju-i{z) + J v+ i(z) = ^J v {z) (42) 



J v -i{z)-J v+l {z) = 2^-J v {z) (43) 
az 



Then 



I — 2^( a ( r / 2 ' Z ) I a(r/2,z)-a(r/2,z'),a(z,z') + a ( r / 2 , z')I a (r /2,z');a(r /2,z) ,a(z,z') + Ol{z, Z > )I a (z,g'y,a(r/2 > z) l a(r/2,g')) (44) 
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where 



I a -p n = I \ i J n +i{ax)J n {(3x)J n {^x)dx (45) 



and Ifl.c, -v and Iy; a ,p are denned as cyclic permutations of the same integral. We use the formula of Sonine and 
Dougall [21| to calculate I a ;p,T It says that for any positive constants a, b and c 



Jn(at) J v (bt) J v (cf )i 1_/x <# 



aPTbi - u)T(y + l/2)r(l/2) J 



(a 2 -b 2 -c 2 + 26c cos p)" - " -1 sin 2 " ip dip (46) 



where 



if a 2 < (b - c) 2 

A a; b,c = { arccos h -^^ if {b - c) 2 < a 2 <(b + c) 2 (47) 

\2 ^ „2 



Thus 



if (b + c) 2 < a 2 



I a( r/2, z ) M r/2, z >)Mz,*') = a{r /2, z)^T(n + l/2)T(l/2) J Sm f d(fi (48) 



We have thus reduced x( r ) to a two dimensional integral: 

2tt 2 ™ 



nT(n+ l/2)r(n)r(l/2) 



l-r/2 /-l-r/2 Mo(f/s,,);«(r/2,,'), («,,') 

cfe / dz'(a(r/2,z')/27r) 2,l (a(z,z')/27r) 2 " / dpsin 2n <p 



l-r/2 /-l-r/2 A(.,<')i«(r/2,!),«Wi,.') 

cfe / dz'(a(r/2,z)/2vr) 2n (a(r/2,z')/27r) 2 ™ / dv?sin 2l > 

II Jo Jo 

l-r/2 /-l-r/2 A(r/l,«),a(r/S l »'),«(',-«') 

cfe / dz'(a(r/2,z')/27r) 2 ™(a(z,-z')/27r) 2 ™ / d^sin 2 ™^ 

o Jo Jo 

l-r/2 /-l-r/2 /•A»(»,-»')i<»Cr/2,*),a(r/2,,') 

cfe / dz'(a(r/2,z)/27r) 2n (a(r/2,z')/27r) 2n / d^sin 2 "^) (49) 

Jo Jo 

(The integral over (/? may be evaluated using We need to determine which values of z and z' correspond to 

which functional form of A a -^ yC . We will use the fact that for all z, z' for which < z, z' < 1 — r/2 

a(r/2, zf < (a(r/2, z') + a(z, z')) 2 (50) 

and 

a{z, z') 2 > {air/2, z) - a(r/2, z 1 )) 2 (51) 

Since z' < 1 — r/2 < r/2, the first inequality is obvious. The second inequality follows from the first inequality. Since 
7(2, —z') > 7(2, z') for all z and z', 7(2, z') could be replaced by 7(2, — z') in 15U|) and (|51|) . It follows from l|47|l . 15U|) 
and ijSTJ) that 



A«rWrtrMrt«*) - I arccos "(-/ 2 ^;);+^;ff-) r/2 ' Z)2 if («W 2 . z ') ^ ^ z ')f < °( r A ^ 

[ arccos H(rA|)^^ 

A a(z , z0;a(r / 2lZ ), a(r/2 , z0 - !f f^/O ^ ^/o / , ,M2 (53) 
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We need to translate the equation a(r/2,z) 2 — (a(r/2,z') — a(z,z')) into an equation involving z and z' . This 
can be done by using the definition of a and expanding both sides. In this way it can be shown that 



if z' > a r (z) 

arccos 



A*(r/2,z);a(r/2,2'),«(z,z') ~ \ „„„„„„ a(r /2,z') 2 +q(z,z' ') 2 ~a(r /2 : z) 2 .,-./. ,., ( " >4 1 

and 



■a(z,z , -a(r/z,z) -r i / x 

2a(r/2,z>)a(z,z>) 112 ^ u r\Z) 



, _ i ■■ if > a r(2) 

\( V ')i«W^),«W2,2') - i a (r/2,z) 2 +a(r/2,z>) 2 -a( Z ,z') 2 . f j ^ „ ,^ 

2a(r 



arccos — ■ — ■ — tt~? — 77; — \ — t — u\ — t\ if z <C & r (z*\ 

°~ (r /2,z)a(r /2,z') ' ^ ' 



where z' = a r (z) is the positive root of the equation 

3 - r 2 - 4z 2 - 4zz' - 4z' 2 - 2rz' + 4r 2 zz' + 8rz V + 8rzz' 2 - 2rz = (56) 

When r = 1 this equation can be factorised as 

(l-2z){l-2z')(l+z + z')=0 (57) 

Hence z' is undetermined whenever z = 1/2 in this case. 
It can be shown in the same way that 

Jo if z' > b r (z) 

A a (r/2,z); a (r/2,z>), a (z,-z') ~ \ ar . nn „ a a(r j2,z') 2 ,- z' ) 2 -q(r /2,z) 2 . f , , , (58) 



arccos " v ' '"" ' , ?} z ), V — a ' Z if z ' < &r(z) 



and 



if z' > b r (z) 

A a ( z ,- z >) M r/2,z)Mr/2.z>) = \ a(r /2,z) 2 + a (r /2,z') 2 - a (z,-z') 2 if v , ^ h ( 59 ) 

2a{r 



arCCOS «~(r/2,z)a(r/2,z') HZ<O r (Z) 



where z' — b r {z) is the positive root of the equation 

2rz + 2rz - 4zz' - 3 + 4z 2 + 4z' 2 + r 2 = (60) 
Since (|56|l and (|60l) are both symmetric in z and z', we could equally well write their solutions as z = a(z') and 



z 



= b(z') instead. a r (z) and b r (z) for r > 1 are shown in figure ^ 
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1/2 



cii(z) 




1/2 



FIG. 2: The functions oi and bi (here r = 1). &i(z) = -! + § + | - 2z)(3 + 2s). 61(0) = when z = ± 



So 



2vr 2 



/■arccos(y a(r/2 , (r/2 , 
\2n / 1, _ • 2n 



nr(n + l/2)r(n)r(l/2) 

/■l-r/2 ,a r (.0 
2 / dz / dz'(a(r/2,z')/27r) 2n (a(z,z , )/27r) 

Jo Jo 

l-r/2 /•a r (*0 ^arccos(y a(i , a >) ;a ( r / 2 ,*),a(r/2 ,»')) 

dz / dz'(a(r/2,2)/27r) 2 "(a(r/2,z')/27r) 2n / d^ sin 2 " 

Jo Jo 

/■l-r/2 rb r (z) r a - 1 ' CCOS (y a (T/2,z);c l ( r /2,z>), a (z,-z>)) 

2 dz dz'(a{r/2,z')/2Tr) 2n (a(z,-z')/2Tt) 2n 

Jo Jo Jo 

l-r/2 fb r (z) (Va(»,-«0!<»Cr/3,«),a(r/a,»0) 

dz / dz'(a(r/2,z)/27r) 2n (a(r/2,z')/27r) 2n / d^ sin 2 " <p 

n Jo JO 

l-r/2 /.1/2 /-tt 

dz / dz'(a(r/2,z)/27r) 2 "(a(r/2,z')/27r) 2 " / sin 2 " ^) 

n J6 r (z) JO 



dp sin " 



7 • 2n 

a</5 sm (p 



(61) 



where y a;/3 
eEl ] 

X(l) 



- + 2p 7 a ■ As r tends to 1, it follows from (|57|l that a r (z) takes the value 1/2 for all z, as shown in 



figure [3 In the special case r — 1 the integral simplifies to 

2?r 2 " 



/•arccos (2/0,(1/2, *) i0 <(l/2,«'), <*(>,*')) 
^2n / J,„ „-2n . 



dip sin 99 



nr(n + l/2)r(n)r(l/2) 

r l/2 rl/2 

x( 2 / dz dz'(a(l/2,z')/27r) 2 "(a(z,z')/27r) s 
Jo Jo 

Z" 1 / 2 /-1/2 / " arCCOS (a Q ( i , i ');c«(l/2. i ). Q (l/2, i ')) 

+ / dz / dz'(a(l/2,z)/27r) 2 "(a(l/2,z')/27r) 2 " / dp sin 2 " ^ 

Jo Jo JO 

rl/2 r&i(z) rarccos (y o(1/2 ,,) i0 (i/2 ,»').«(* ,-*')) 

+ 2 / dz dz'(a(l/2,z')/27r) 2 "(a(z,-z')/27r) 2 " / d^ sin 2 " p 

Jo Jo Jo 

1/2 fbl(z) r a - ICC ° S (V a (z,-z'); a (l/2,z), a (l/2,z')) 

dz / dz'(a(l/2,z)/27r) 2 "(a(l/2,z')/27r) 2 " / dp sin 2 " ^ 

• Jo Jo 

1/2 ,1/2 /.7T 

-!- / dz / dz'(a(l/2,z)/27r) 2 "(a(l/2,z')/27r) 2 " / d^ sin 2 " p) (62) 

Jbi(z) JO 



After integration by parts, this gives integrals of the type 

p(x) 



I 



q{x)\/ a + bx + ex 1 



-dx 



(63) 



10 



where p and q are polynomials. Using Maple it was thus possible to calculate f° r D = 5, 7, 9, 11. We may 
now obtain [/^ from 1)17(1 . Since Q and [/] can be obtained from l(20|l and 121(1 . we have found B4. We use the 
more compact Ree Hoover / formalism j23] to present the results. In this formalism B4 consists of only two diagrams 
instead of three. Here 

/(r«) - /(r«) = 1 (64) 

Thus 

= □ 

! = S + 2D+D (65) 

and 

^ 4 = I ~ I 1 (66) 

The final answer is given in tables ITT1 IIIII and IIVI The numerical values of references 0] and 0] agree with the exact 
result. 



TABLE III: Exact and numerical J, 7] values of the Ree Hoover complete star 



D 





decimal expansion 


3 


89 219 s/2 | 4131 arccos (1/3) 
280 1120 tt ' 2240 n 


0.31672598803... 
0.31673(2) 


5 


163547 3888425 V2 , 67183425 arccos (1/3) 
128128 8200192 tt 1 10400384 tt 


0.11520591833... 
0.115211(3) 


7 


283003297 159966456685 \/2 , 292926667005 arccos (1/3) 
141892608 217947045888 tt 1 48432676864 tt 


0.04492254969... 
0.044927(2) 


9 


88041062201 2698457589952103 V2 . 8656066770083523 arccos (1/3) 
34810986496 2851716013752320 tt 1 1140686405500928 rr 


0.01828214224... 
0.018286(1) 


11 


66555106087399 16554115383300832799 s/2 , 52251492946866520923 arccos (1/3) 


0.00766164876... 
0.0076638(8) 


22760055898112 14916030233386680320 tt 1 5966412093354672128 tt 



TABLE IV: Exact and numerical 17) values of the Ree Hoover ring 



D 
3 

5 

7 

9 

11 



3 

' SB'i 



decimal expansion 



4131 , 657 V2 _ 12393 arccos (1/3) 
4480 ' 2240 tt 4480 tt 



67183425 , 11665275 s/2 _ 201550275 arccos (1/3) 
32800768 ' 16400384 n 32800768 tt 



292926667005 , ir,'i')i;[,-l.M;(,8r, ^2 _ 878780001015 arccos (1/3) 
9(;.s(;r,35:i728 + 145298030592 tt 96865353728 tt 



8656066770083523 , 8095372769856309 V2 _ 25968200310250569 arccos (1/3) 
2281372811001856 5703432027504640 tt 2281372811001850 tt 



52251492946866520923 , 49662346149902498397 %/2 _ 156754478840599562769 arccos (1/3) 
11932824186709344256 ' 29832060466773360640 tt 11932824180709344256 tt 



-0.02977648205.. 

-0.029781(8) 

-0.03923343804.. 

-0.039233(3) 

-0.03505760307.. 

-0.035055(3) 

-0.02686294042.. 

-0.026861(3) 

-0.01899884734.. 

-0.018997(3) 



IV. DISCUSSION 



Table [H] shows exact and numerical values of the fourth virial coefficient. Typically the relative error of the 
numerical value is of order 10 -4 . Recently Clisby and McCoy 01 calculated higher order coefficients using Monte 
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TABLE V: Contributions to the fifth virial coefficient in three dimensions 123 



diagram 



exact value 



E5/B% 

E&a/Bt 

E6f3/B% 

E7a/B% 

E7f3/Bi 

E7 1 /B$ 

E8a/B% 

E8p/B$ 

E9/BI 

EW/Bl 



10752 

68419 

26SS0 

82 

35 

34133 
" 17920 

18583 
" 5376 

73491 



+ 



33291 y/5" 
9800 7T 



unknown 

35731 i 1458339 V2 33291 V3 , 683559 arccos (1/3) 



6720 627200 

unknown 
unknown 



9800 



.'',5840 



Carlo methods. It is seen that the relative error increases with the order of the coefficient, which is one of the 
reasons why it is desirable to find analytic values. Table shows the known exact values of diagrams of the fifth 
virial coefficient in 3 dimensions. The diagrams E7(3 and E8f3 have the same coefficient, so there is so far no total 
contribution of %/3/tt- We make the following conjecture: 

Conjecture 1 In any dimension, the hard sphere potential allows the analytic computation of every virial coef- 
ficient B n . 

No one has so far been able to prove this conjecture, but since the complete star of any number of points can be 
expressed in terms of an integral involving the same functions F and G used in section III, there is good reason to 
assume that it is true and can be proven. 
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